Rectangular quantum dots in high magnetic fields 
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We use density-functional methods to study the effects of an external magnetic field on two- 
dimensional quantum dots with a rectangular hard-wall confining potential. The increasing magnetic 
field leads to spin polarization and formation of a highly inhomogeneous maximum-density droplet 
at the predicted magnetic field strength. At higher fields, we find an oscillating behavior in the 
electron density and in the magnetization of the dot. We identify a rich variety of phenomena 
behind the periodicity and analyze the complicated many-electron dynamics, which is shown to be 
highly dependent on the shape of the quantum dot. 

PACS numbers: 73.21. La, 71.10.-w 



I. INTRODUCTION 

The problem of many interacting electrons in a non- 
circular quantum well under the influence of an external 
magnetic field is extremely challenging for the compu- 
tational tools available. For the time being, theoretical 
studies of such systems have dealt either with very small 
electron numbers, or with the single-electron properties 
in non-integrable quantum billiardsi On the other hand, 
the rapid technical development in the realization of dif- 
ferent quantum-dot systems naturally motivates theoret- 
ical modeling of complicated many-electron structures.^ 

The magnetic-field dependence of the energy spectrum 
in a square two-electron quantum dot has been analyzed 
in detail by Creffield et al^ They found Aharonov-Bohm- 
type oscillations in the lowest levels, indicating periodic 
singlet-triplet changes in the ground state. Ugajini found 
that these state transitions lead to strong effects in the 
optical excitation spectra. Within a similar square-dot 
system, he studied the effects of the Coulomb interaction 
on the far-infrared-absorption (FIR) spectra*^ Recent 
density-functional FIR calculations for soft- wall triangu- 
lar and square dots have been done by Vali'n-Rodriquez 
et al.f' who identified corner and side modes in the sys- 
tem. The interactions and the quantum-dot geometry 
affect also strongly the magnetization, at least for small 
electron numbers 

In a circular geometry, a maximum-density droplet 
(MDD) caused by the magnetic field corresponds to a 
polarized state with electrons occupying successive an- 
gular momentum states, i.e., I — 0, — 1, —N + 1, giv- 
ing Lmdd — —\N{N — 1) for the total angular momen- 
tum of N electrons. A stability region for the MDD 
has been identified in electron transport experiments 
through vertical quantum dots by Oosterkamp et al& 
This state, as well as the post-MDD regime in higher 
magnetic fields, has been considered theoretically in cir- 
cular parabolically-confined quantum dots by several au- 
|-j-^Qj,g^i0jiiji2 ^Yie case of a non-circular symmetry, how- 
ever, the high-magnetic-field limit has not been previ- 
ously analyzed to the best of our knowledge. 

In this work, we examine the magnetic-field effects 
in rectangular quantum dots by using the spin-density- 



functional theory (SDFT). Our aim is to define the MDD 
state in a rectangular geometry and extend our previous 
analysis of the MDD formation in hard-wall quantum 
dots.^^ The onset of the MDD can be predicted from the 
number of flux quanta which only depends on the area of 
the dot and is thus geometry- independent. We study also 
the beyond-MDD regime that can be characterized by pe- 
riodic state oscillations. The periodicity can be observed 
in the magnetization which we compare with the cor- 
responding results for noninteracting electron systems. 
The origin of the periodicity is generally the competition 
between magnetic confinement and Coulomb repulsion, 
but the behavior of the single-electron states is very sen- 
sitive to the dot geometry. 

This paper is organized as follows. In Sec^Jwe present 
briefly the model Hamiltonian and the computational 
methods, based on a real-space spin-density-functional 
approach. In Sec. IIIII we consider the formation and the 
structure of the MDD in rectangular geometries. Sec- 
tion IIVI presents the state oscillations above the MDD 
regime, first from the point of view of the total magne- 
tization, and finally at the level of the effective single- 
electron states. The paper is summarized in Sec. 



II. THE MODEL AND THE METHODS 

We restrict the dot to the xy plane and use the effec- 
tive mass approximation with the material parameters 
for GaAs, i.e., the effective mass m*=0.067 me and the 
dielectric constant e = 12.7. The many-body Hamilto- 
nian reads 



H = 



1 ^ 
2m* ^ I 



N 



nVi + -A(r,) + ^ 
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(1) 
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where the shape of the dot is defined by a rectangular 
hard-wall confining potential. 



0, 0<x<(3L,0<y<L , , 
oo, otherwise. ^ ' 
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Here [3 defines the ratio of the side legths of the rectangle, 
and L is scaled such that the area of the dot remains con- 
stant, i.e., I3L'^ = 47r^a^^, where la*gK. 10.03 nm. We 
apply an external magnetic field B — Bz perpendicular 
on the xy plane, and use the symmetric gauge, 



A = — (-y,x,0), 



(3) 



for the vector potential. The last term in Eq. 0J is the 
Zeeman energy (with g* — —0.44 for GaAs) which has 
an insignificant effect on the results presented and is thus 
neglected. In the calculations we employ the SDFT in 
the self-consistent Kohn-Sham (KS) formulation. In the 
systems considered in this study, the SDFT gives rea- 
sonably accurate results compared to the computation- 
ally more demanding current-spin-density-functional ap- 
proach. In the local spin-density approximation we use 
the exchange-correlation energy by Attaccalite et alm^ 
It is based on the diffusion Monte Carlo simulations over 
the whole range of spin polarization, which is a major im- 
provement comparing to the previous parametrizations. 
This leads to more accurate results in both zero and non- 
zero magnetic fields as shown in Ref. fl^ 

The calculations are perfomed in real space using fi- 
nite differences for the derivative operations on two- 
dimensional point grids. Since there are no implicit re- 
strictions for the symmetry, the external potential can be 
shaped arbitrarily in the computing region. The number 
of grid points is 128 x 128, giving a sufficient accuracy 
in the total energy. To make the convergence effective, 
we employ the Rayleigh quotient multigrid methociiS for 
the discretized single-electron Schrodinger equations. 



III. MDD FORMATION 
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FIG. 1: Chemical potentials of the fully polarized states for 
A^-electron rectangular quantum dots as a function of the 
magnetic field B. 
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FIG. 2: MDD-window limits as a function of the number of 
electrons in different quantum-dot geometries. The line for 
the predicted MDD formation is also shown. 



In the zero-field solution for a rectangular quantum dot 
with an even number of electrons, the total spin is cither 
zero or one, depending on the degeneracy on the highest 
leveliii Increasing B brings the higher spin states closer 
to the S* = and = 1 states in energy, and finally the 
full spin polarization is found. The value of B for the 
complete spin polarization is rather independent of the 
deformation parameterii^ In the eight-electron dot, for 
example, the ground state becomes fully spin-polarized 
at B - 10 T as = 1 and at 11 T as /3 = 2. 

Figure ^ shows the chemical potentials, n{N) — 
E{N) — E{N — 1), of the polarized states as a function of 
B for TV = 7 — 12. The point when the totally polarized 
state becomes the ground state is marked with a dotted 
line in the figure. We consider only the polarized state in 
order to clarify the MDD window, which corresponds to 
the descending regime in the chemical potential. Similar 
behavior in /i was observed in the experiments for vertical 
quantum dots by by Oosterkamp et al^. They measured 
the evolution of the Coulomb blockade peaks as a func- 
tion of S for = — 40 and found clearly identifiable 
phases for the filling factor, of which v = 1 corresponds 



to a maximum-density droplet. According to our defi- 
nition, the true MDD exists on the rightmost ascending 
stripe in Fig. ^ As /3 increases, the MDD windows be- 
come larger and flatter due to pronounced localization 
in the x direction. Thus, the square case shows the most 
drastic behavior because of the highest consistency of the 
dot symmetry with the magnetic confinement. 

Figure [21 gives the limits for the MDD window as a 
function of N in different quantum-dot geometries. We 
see that the increase in the width of the MDD window 
is due to a shift in the reconstruction point, whereas the 
onset of the window is nearly independent of the dot 
geometry. The reason is the following. The formation of 
the MDD requires an equal number of fiux quanta, N,^, — 
<I>/<I>o, to the A^ — 1 vortices 'seen' by an electron as the 
other electrons are fixed. A vortex is defined as a point 
where the wave function is zero and its phase changes by 
—2tt for a counterclockwise rotationjl^ Here $o = h/e 
is the single flux quantum, and since the magnetic flux 
is directly proportional to the dot area, i.e., <I> = BA^ 
we thus get B = h{N - l)/eA « 1.04(A^ - 1) T for 
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FIG. 3: Electron densities for the MDD states of eight- 
electron quantum dots with different geometries. 



the onset of the MDD in all geometries. We show this 
prediction in Fig.[21and it agrees remarkably well with the 
transition points obtained from the chemical potentials. 
The reconstruction limit is, however, much more difficult 
to estimate and it depends strongly on the geometry as 
explained above. 

At the end of the MDD window the droplet achieves 
the smoothest electron density. In Fig. |31 we show the 
maximum- density distributions for /3 = 1 and /3 = 2 
quantum dots, and also for a circular hard-wall dot with 
N = 8. They have a highly inhomogeneous structure 
instead of a flat distribution found in parabolically con- 
fined MDD'siii This is a direct consequence of the hard- 
wall confinement where the Coulomb interaction favors 
strong localization in the corners and on the edges over 
the whole range of B. In the circular hard-wall MDD, 
the occupation of the angular momentum orbitals follows 
the original definition of a maximum-density droplet for 
a parabolic dot, i.e., Iz = 0,-1,...,— 7. The number of 
vortices in a particular KS state thus equals \lz\. For N 
electrons at the positions Zi — Xi + iyi, the many-body 
wave function for a parabolic MDD can be written in 
Laughlin's^^ form as 

^MDD = exp ( -|b$o J2 i^'^n " ^j)' ("^^ 

V k / i<j 

which is similar to the wave function for a state with the 
filling factor v — I in the thermodynamic limit. It is easy 
to see that by following the above definition for a vortex 
we obtain the correct total number of iV — 1 vortices in 
the system. 

The rectangular dot shows a more complicated struc- 
ture: the formation of the MDD corresponds to the ap- 
pearance of the correct maximum number of iV — 1 vor- 
tices in the KS state highest in energy, but in the other 
states there are typically 'extra' vortices, leading to more 
than N{N — l)/2 of them in the KS states altogether. 
The reason is the complex structure of the indefinite an- 
gular momentum states that prevents the description of 
the non-circular MDD with a Laughlin-type many-body 
wave function, even if the B-N dependence for the MDD 
formation can be correctly estimated. Thus, we resort 
to considering the evolution of the expectation values of 
the angular momentum operator, Iz = —ih{x-^ — y^), 
for different KS states in the high-i? regime. Below we 
designate this expectation value for non-circular dots as 
the effective Iz- 



IV. STATE OSCILLATIONS 

The rapidly growing Coulomb interaction leads to the 
reconstruction of the MDD into a lower-density droplet 
with a higher angular momentum as B increases. Ac- 
cording to our calculations, the collapse of the circu- 
lar hard-wall MDD leads into a ringlikc structure with 
Iz = — 1, — 2, — iV, contrary to the mean-field results 
for parabolic dots, where the Chamon-Wen edge states 
and charge-density waves have been foundi^i Due to the 
hard- wall confinement in our case, there is a remarkably 
higher density on the edge already in the MDD state, 
and the exchange-correlation energy is thus relatively less 
important than in a parabolic dot. However, exact di- 
agonalization and quantum Monte Carlo results predict 
circularly symmetric reconstruction also in the parabolic 
case^ In a rectangular geometry, the result is a strongly 
localized edge state that resembles a Wigner molecule in 
the sense that there are pronounced peaks in the electron 
density to minimize the Coulomb interaction^ 

As B is increased further, the electronic structure be- 
gins to oscillate between well-localized and more diffuse 
states. This is a result from the interplay between the 
magnetic confinement and the Coulomb interaction. In 
this connection, the polarized ground state changes peri- 
odically as the interactions entangle the single-electron 
states on the lowest Landau level. This is similar 
to singlet-triplet oscillations found in interacting two- 
electron quantum dots of both circular^^ and square^ 
shapes, as well as in a double-dot systemiS^ 



A. Magnetization 

The structural oscillations beyond the MDD state can 
be observed in the magnetization of the dot, since the 
growth in the total angular momentum leads to the de- 
crease in the current. The total magnetization for a finite 
system is defined at zero temperature as 



In Fig. 21 we show the magnetization in three rectangu- 
lar quantum dots beyond the MDD regime. The local 
minima correspond to the points where the localization 
on the edge is maximized, and at local maxima there are 
states with a density contribution at the center. It can 
be clearly seen in Fig.0]that M decreases with increasing 
electron number, and the onset of the oscillation moves to 
higher B as the deformation f3 is increased, which is con- 
sistent with the upper limit of the MDD window shown 
in Fig. 121 The period in magnetic field strength is nearly 
a constant, ^ 2.5 T in all the systems presented in Fig.^ 
We emphasize that the majority of (/?, N) combinations 
lead to more complicated and even irregular magnetiza- 
tion. This will be discussed in more detail in Sec. HVBl 
Figure^ also shows the noninteracting magnetization for 




{(3,N) — (2,6), calculated from the six lowest eigenener- 
gies on the first Landau level. In the MDD regime, it is 

1 fig above the interacting curve with the same slope. 
At i? ~ 15 T, however, there is a clear crossing which 
indicates that the interactions retard the overall increase 
in \M\ as a function on B in the oscillatory regime. 

Magnetization in square-shaped quantum dots up to 
iV = 4 has been studied by Magniisdottir and Gud- 
mundssoni They found remarkably higher magnetiza- 
tion in the noninteracting than in the interacting pic- 
ture, but did not consider the high-_B limit, where we see 
the crossing between those two cases. In a two-electron 
square dot studied by Sheng and Xu,^ the magnetization 
looks similar to our case, but the oscillation represents 
the trivial singlet-triplet exchange4i 

B. Development of the angular momenta 

As a mean-field theory, the density-functional ap- 
proach only approximates the many-body energy states. 
However, we will next consider the development of the 
effective angular momentum of different KS orbitals in 
order to obtain a deeper insight on the state oscillations 
on the single-electron level. In Fig. [S] we present the evo- 
lution of the effective Iz values of the occupied KS levels 
in the polarized state of a square eight-electron quantum 
dot as a function of B. The structure is intriguing and 
shows clearly two different modes in the system. First, 
there are four closely entangled states, representing the 
relatively stable corner modes. They are reminiscent of 
the Aharonov-Bohm-like energy-level structure in an in- 
teracting two-electron dot,- Second, there are states that 
with increasing B condense to nearly circular orbits on 
the first Landau level before the Coulomb interaction 
expands them to high effective Iz values. In this con- 
text, the state development is similar to the circular case. 



FIG. 5: Evolution of the effective Iz values as a function of 
the magnetic field B in a square eight-electron quantum dot. 

where the occupation shifts from the lowest Iz state to 
the highest level as the ground state changes. The pe- 
riodic spread of the Landau states is clearly the origin 
of the magnetization oscillations discussed above. The 
upper limit of the MDD window, marked in Fig. [S] can 
be interpreted as the most evenly extended effective Iz 
distribution similar to the corresponding circularly sym- 
metric system. 

Figure ini shows the development of the effective Iz val- 
ues in the case oi f3 — 2. Due to the reduction of the 
symmetry, the evolution is more complicated than in 
the square case, and we can find four different modes. 
First, there are two corner modes in the system. Of 
them, the highest effective Iz state shows steady local- 
ization at the both ends of the rectangle, and the second 
one (5) oscillates weakly. The latter is visualized in the 
right panel of Fig . El with the KS wave functions, IV'ksI^ 
(see also Ref. l23l) . In addition, there are two different 
types of highly oscillatory states. The first pair (3,4) has 
a remarkably lower amplitude in the oscillation regime 
{B > 14 T). As can be seen in the KS wave functions, 
they correspond to localization in the y direction. Inter- 
locked to these states, there are strongly oscillative states 
(1,2) that, according to the KS wave functions, seem to 
form a periodic structure in the rectangle, reminiscent of 
a 'lozenge' orbit in the correponding classical bouncing 
map, studied in the square case by AguiariSi It is notice- 
able that these states affect strongly on the corner modes 
by compressing them slightly to lower Iz- The cusps at 
B = 11.8 T and 13.8 T are numerical noise, correspond- 
ing to degeneracies of the KS orbitals 1 and 3 or 2 and 
4, respectively. 

It is clear from Fig. El that the formation of bulk Lan- 
dau states disappears as the dot is squeezed toward a 
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FIG. 6: Evolution of the effective values as a function of the magnetic field B in a P — 2 eight-electron quantum dot. KS 
wave functions, corresponding to different modes, are also shown. 



rectangular shape. Thus, the electron-electron interac- 
tions make the quantum dot sensitive to the deformation 
parameter. On the other hand, the overall behavior be- 
yond the MDD is qualitatively similar in all geometries, 
even if the evolution of the effective Iz values gets very 
complicated. This occurs e.g. in the case of /3 = 1.5, 
where the different modes discussed above are hardly rec- 
ognizable. In this sense, an increase in the symmetry sep- 
arates the different states and makes the corresponding 
classically stable orbits visible. 

We underline that within the density-functional the- 
ory, one always has to regard carefully the KS states and 
the corresponding wave functions shown in Figs. |S1 and 
In this particular system, however, we find that the pro- 
cedure gives important information on the many-particle 
dynamics in an interacting non-intcgrablc quantum dot, 
not accessible with more 'exact' computational methods. 

V. SUMMARY 

To summarize, we have studied the electronic struc- 
ture of rectangular quantum dots in a magnetic field. 
We have identified the formation of a maximum-density 
droplet at the predicted B, even if the electron density is 
strongly localized due to the hard- wall confinement. Be- 



yond the MDD, we have found oscillating magnetization 
that reflects the periodic changes in the polarized ground 
state. In the oscillating regime, the magnetization shifts 
to lower values than in the corresponding noninteracting 
case. The origin of the periodic behavior has been ana- 
lyzed by following the evolution of the effective angular 
momentum quantum numbers. There are clear corner 
modes in all rectangular systems, but the development 
of the center orbits is found to be very sensitive to the 
geometry of the dot. In a square, the central KS orbitals 
are nearly symmetric, whereas a rectangle with (3 = 2 
shows two different oscillative modes. When the symme- 
try is reduced further, the center orbitals become unin- 
dentifiable. In spite of the limitations in the mean-field 
approach, we have managed to bring out the richness of 
phenomena exhibited by the electronic structure in the 
high-i? limit. We hope that this would inspire new ex- 
perimental work on this topic. 
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